A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.

Intercontinental dispersal prior to human translocation revealed in a cryptogenic invasive tree

In this study, complementary species-level and intraspecific phylogenies were used to better circumscribe the original native range and history of translocation of the invasive tree Parkinsonia aculeata. Species-level phylogenies were reconstructed using three chloroplast gene regions, and amplified fragment length polymorphism (AFLP) markers were used to reconstruct the intraspecific phylogeny. Together, these phylogenies revealed the timescale of transcontinental lineage divergence and the likely source of recent introductions of the invasive. The sequence data showed that divergence between North American and Argentinean P. aculeata occurred at least 5.7 million years ago, refuting previous hypotheses of recent dispersal between North and South America. AFLP phylogenies revealed the most likely sources of naturalized populations. The AFLP data also identified putatively introgressed plants, underlining the importance of wide sampling of AFLPs and of comparison with uniparentally inherited marker data when investigating hybridizing groups. Although P. aculeata has generally been considered North American, these data show that the original native range of P. aculeata included South America; recent introductions to Africa and Australia are most likely to have occurred from South American populations.

Is a native rodent competitively dominant over an invasive rodent in lowland agro-forest habitat of the Philippines?

In the lowland agro-forest of the Sierra Madre Biodiversity Corridor (SMBC), it is considered that a native rodent species, Rattus everetti is competitively dominant over an invasive pest species, Rattus tanezumi. The main aim of this study was to assess the response of R. tanezumi following short term removal of R. everetti. We tested this experimentally by trapping and removing R. everetti from two treatment sites in agro-forest habitat on three occasions over three consecutive months. This was followed by three months of non-removal trapping. Two non-treatment sites were trapped for comparison. Following R. everetti removal, R. everetti individuals rapidly immigrated into the treatment sites and a significantly higher proportion of R. tanezumi females were in breeding condition in the treatment sites than in the non-treatment sites. The results from this study provide evidence of competition between native and invasive rodent species in complex agro-ecosystems. We were also able to demonstrate that R. everetti populations are able to recover rapidly from the non-target effects of short-term lethal control in and around agro-forest.

Selective toxicity of different granular insecticides to the peach-potato aphid, Myzus persicae (Homoptera: Aphididae) and its parasitoid, Aphidius matricariae (Hymenoptera: Aphidiidae) in two differentially resistant potato cultivars

Field studies were conducted in Pakistan to examine the effects and the interaction of two differentially resistant potato cultivars i.e. Cardinal and Desiree (one partially resistant and one susceptible to Myzus persicae (Sulzer), respectively) with different dosage rates of granular insecticides, at different time intervals after application in relation to percent kill of M. persicae and effects on the parasitoid Aphidius matricariae Haliday (i.e. the third trophic level) within the aphid mummies, percent parasitism and Thimet 10G (phorate) was found about 30% more effective in reducing aphid population than the Furadan 3G (carbofuran). The highest doses of each insecticide caused maximum aphid mortality, in general aphid mortality appeared dose dependent. Mostly all the higher and lower doses of the tested insecticides were about 10% more effective on Cardinal than on Desiree. The most significant result was the synergistic interaction at the lower doses with plant resistance, so that the same level of control was recorded with the second highest dose on Cardinal as with the highest dose on Desiree. Also the same level of control was observed at the lowest dose on Cardinal as with the second last lowest dose on Desiree. Furadan 3G was found least toxic to the A. matricariae in terms of percent parasitism, emergence of parasitoids and number of mature eggs in the emerging females. Furadan 3G gave 13, 15 and 6% higher figures, respectively from the parasitoid characteristics than Thimet 10G. The highest doses of both insecticides were clearly toxic to the parasitoid. In general, the effects on the parasitoid were dose dependent. The second highest dose of Thimet 10G, gave the maximum yield

Stocking may increase mitochondrial DNA diversity but fails to halt the decline of endangered Atlantic salmon populations

Over the last 50 years, Spanish Atlantic salmon (Salmo salar) populations have been in decline. In order to bolster these populations, rivers were stocked with fish of northern European origin during the period 1974-1996, probably also introducing the furunculosis-inducing pathogen, Aeromonas salmonicida. Here we assess the relative importance of processes influencing mitochondrial (mt)DNA variability in these populations from 1948 to 2002. Genetic material collected over this period from four rivers in northern Spain (Cantabria) was used to detect variability at the mtDNA ND1 gene. Before stocking, a single haplotype was found at high frequency (0.980). Following stocking, haplotype diversity (h) increased in all rivers (mean h before stocking was 0.041, and 0.245 afterwards). These increases were due principally to the dramatic increase in frequency of a previously very low frequency haplotype, reported at higher frequencies in northern European populations proximate to those used to stock Cantabrian rivers. Genetic structuring increased after stocking: among-river differentiation was low before stocking (1950s/1960s Phi(ST) = -0.00296-0.00284), increasing considerably at the height of stocking (1980s Phi(ST) = 0.18932) and decreasing post-stocking (1990s/2002 Phi(ST) = 0.04934-0.03852). Gene flow from stocked fish therefore seems to have had a substantial role in increasing mtDNA variability. Additionally, we found significant differentiation between individuals that had probably died from infectious disease and apparently healthy, angled fish, suggesting a possible role for pathogen-driven selection of mtDNA variation. Our results suggest that stocking with non-native fish may increase genetic diversity in the short term, but may not reverse population declines.

A note on the accuracy of PAC-likelihood inference with microsatellite data

Stephens and Donnelly have introduced a simple yet powerful importance sampling scheme for computing the likelihood in population genetic models. Fundamental to the method is an approximation to the conditional probability of the allelic type of an additional gene, given those currently in the sample. As noted by Li and Stephens, the product of these conditional probabilities for a sequence of draws that gives the frequency of allelic types in a sample is an approximation to the likelihood, and can be used directly in inference. The aim of this note is to demonstrate the high level of accuracy of "product of approximate conditionals" (PAC) likelihood when used with microsatellite data. Results obtained on simulated microsatellite data show that this strategy leads to a negligible bias over a wide range of the scaled mutation parameter theta. Furthermore, the sampling variance of likelihood estimates as well as the computation time are lower than that obtained with importance sampling on the whole range of theta. It follows that this approach represents an efficient substitute to IS algorithms in computer intensive (e.g. MCMC) inference methods in population genetics. (c) 2006 Elsevier Inc. All rights reserved.

Population genetic models can be used to study the evolution of the interacting behaviours of parents and their progeny

An example of the evolution of the interacting behaviours of parents and progeny is studied using iterative equations linking the frequencies of the gametes produced by the progeny to the frequencies of the gametes in the parental generation. This population genetics approach shows that a model in which both behaviours are determined by a single locus can lead to a stable equilibrium in which the two behaviours continue to segregate. A model in which the behaviours are determined by genes at two separate loci leads eventually to fixation of the alleles at both loci but this can take many generations of selection. Models of the type described in this paper will be needed to understand the evolution of complex behaviour when genomic or experimental information is available about the genetic determinants of behaviour and the selective values of different genomes. (c) 2007 Elsevier Inc. All rights reserved.

Contrasting effects of necrotrophic and biotrophic plant pathogens on the aphid Aphis fabae

Phytophagous insects have to contend with a wide variation in food quality brought about by a variety of factors intrinsic and extrinsic to the plant. One of the most important factors is infection by plant pathogenic fungi. Necrotrophic and biotrophic plant pathogenic fungi may have contrasting effects on insect herbivores due to their different infection mechanisms and induction of different resistance pathways, although this has been little studied and there has been no study of their combined effect. We studied the effect of the biotrophic rust fungus Uromyces viciae-fabae (Pers.) Schroet (Basidiomycota: Uredinales: Pucciniaceae) and the necrotrophic fungus Botrytis cinerea Pers. (Ascomycota: Helotiales: Sclerotiniaceae) singly and together on the performance of the aphid Aphis fabae Scop. (Hemiptera: Aphididae) on Vicia faba (L.) (Fabaceae). Alone, botrytis had an inhibitory effect on individual A. fabae development, survival and fecundity, while rust infection consistently enhanced individual aphids’ performance. These effects varied in linear relation to lesion or pustule density. However, whole-plant infection by either pathogen resulted in a smaller aphid population of smaller aphids than on uninfected plants, indicating a lowering of aphid carrying capacity with infection. When both fungi were applied simultaneously to a leaf they generally cancelled the effect of each other out, resulting in most performance parameters being similar to the controls, although fecundity was reduced. However, sequential plant infection (pathogens applied five days apart) led to a 70% decrease in fecundity and 50% reduction in intrinsic rate of increase. The application of rust before botrytis had a greater inhibitory effect on aphids than applying botrytis before rust. Rust infection increased leaf total nitrogen concentration by 30% while infection by botrytis with or without rust led to a 38% decrease. The aphids’ responses to the two plant pathogens individually is consistent with the alteration in plant nutrient content by infection and also the induction of different plant defence pathways and the possible cross-talk between them. This is the first demonstration of the complex effects of the dual infection of a plant by contrasting pathogens on insect herbivores. Key words: Vicia faba, Botrytis cinerea, Uromyces viciae-fabae, tripartite interactions, induced resistance

Within-host evolution of Staphylococcus aureus during asymptomatic carriage

Background Staphylococcus aureus is a major cause of healthcare associated mortality, but like many important bacterial pathogens, it is a common constituent of the normal human body flora. Around a third of healthy adults are carriers. Recent evidence suggests that evolution of S. aureus during nasal carriage may be associated with progression to invasive disease. However, a more detailed understanding of within-host evolution under natural conditions is required to appreciate the evolutionary and mechanistic reasons why commensal bacteria such as S. aureus cause disease. Therefore we examined in detail the evolutionary dynamics of normal, asymptomatic carriage. Sequencing a total of 131 genomes across 13 singly colonized hosts using the Illumina platform, we investigated diversity, selection, population dynamics and transmission during the short-term evolution of S. aureus. Principal Findings We characterized the processes by which the raw material for evolution is generated: micro-mutation (point mutation and small insertions/deletions), macro-mutation (large insertions/deletions) and the loss or acquisition of mobile elements (plasmids and bacteriophages). Through an analysis of synonymous, non-synonymous and intergenic mutations we discovered a fitness landscape dominated by purifying selection, with rare examples of adaptive change in genes encoding surface-anchored proteins and an enterotoxin. We found evidence for dramatic, hundred-fold fluctuations in the size of the within-host population over time, which we related to the cycle of colonization and clearance. Using a newly-developed population genetics approach to detect recent transmission among hosts, we revealed evidence for recent transmission between some of our subjects, including a husband and wife both carrying populations of methicillin-resistant S. aureus (MRSA). Significance This investigation begins to paint a picture of the within-host evolution of an important bacterial pathogen during its prevailing natural state, asymptomatic carriage. These results also have wider significance as a benchmark for future systematic studies of evolution during invasive S. aureus disease.

Is behavioral ecology important for understanding and predicting population dynamics?

Population ecology is a discipline that studies changes in the number and composition (age, sex) of the individuals that form a population. Many of the mechanisms that generate these changes are associated with individual behavior, for example how individuals defend their territories, find mates or disperse. Therefore, it is important to model population dynamics considering the potential influence of behavior on the modeled dynamics. This study illustrates the diversity of behaviors that influence population dynamics describing several methods that allow integrating behavior into population models and range from simpler models that only consider the number of individuals to complex individual-based models that capture great levels of detail. A series of examples shows the importance of explicitly considering behavior in population modeling to avoid reaching erroneous conclusions. This integration is particularly relevant for conservation, as incorrect predictions regarding the dynamics of populations of conservation interest can lead to inadequate assessment and management. Improved predictions can favor effective protection of species and better use of the limited financial and human conservation resources.

Adaptive approximate Bayesian computation

Sequential techniques can enhance the efficiency of the approximate Bayesian computation algorithm, as in Sisson et al.'s (2007) partial rejection control version. While this method is based upon the theoretical works of Del Moral et al. (2006), the application to approximate Bayesian computation results in a bias in the approximation to the posterior. An alternative version based on genuine importance sampling arguments bypasses this difficulty, in connection with the population Monte Carlo method of Cappe et al. (2004), and it includes an automatic scaling of the forward kernel. When applied to a population genetics example, it compares favourably with two other versions of the approximate algorithm.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.

Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Evidence for a residual postglacial founder effect in a highly dispersive freshwater invertebrate

In the northern hemisphere, glacial periods have had profound and lasting effects on the population genetics of numerous species, with founder effects often persisting for many generations in formerly glaciated regions. We found an unusual example of this in the freshwater bryozoan Cristatella mucedo, which showed regional differences in haplotype diversity with relatively low levels of haplotype diversity in northern Europe compared to central/southern Europe despite previous evidence for frequent dispersal between C. mucedo populations. Such contradictions between high dispersal and low gene flow have now been reported in several other freshwater taxa and may be attributed to persistent founder effects following colonization of sites by a few individuals whose efficient reproduction leads to rapid population growth. Alternatively, selection may determine which genotypes can thrive in northerly locations, or it may be that C. mucedo has undergone cryptic speciation. Future work on adaptive genomic regions is required before we can understand how gene flow, local adaptation, and speciation influence the current distribution patterns of bryozoans and other freshwater invertebrates.